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Abstract 



We present a variational approach for directed polymers in D transversal dimensions 
^ ' which is used to compute the corrections to the mean field theory predictions with broken 

O ■ replica symmetry. The trial function is taken to be a symmetrized version of the mean- 

field solution, which is known to be exact for D = oo. We compute the free energy 
corresponding to that function and show that the finite- Z? corrections behave like D~^'^. 
^ ' It means that the expansion in powers of 1/D should be used with great care here. We 

hope that the techniques developed in this note will be useful also in the study of spin 
glasses. 

1 Introduction 

At the present moment a very useful approach in the study of phase transitions in disordered 
systems is given by the mean-field theory in the framework of the spontaneous breaking of 
the replica symmetry |l| . In many cases there is a well-developed theoretical understanding of 
the mean-field theory, using the replica formalism or using the equivalent probabilistic cavity 
approach. The situation is not so good if we consider the perturbative corrections to the mean 
field approach in the phase where replica symmetry is broken: a few computations exist, but 



they are very difficult from the technical point of view, so that only one loop contributions 
have been computed [H- 

The aim of this paper is to explore a new field: non perturbative effects. Only very few 
papers exist which are devoted to the study of non perturbative effects, among them we recall 
B S 0- ^°^ perturbative effects are very important because they could spoil the predictions 
of mean field theory in a rather subtle way, which cannot be detected by the usual perturbative 
methods. 

In this note we will concentrate our attention on the following problem. In the phase 
where replica symmetry is broken, in replica space there are many different stationary points 
for the free energy as a function of the order parameter (let us call it Qa^b and denote its 
average in a given stationary point by qa^h) ■ These stationary points are related to the other by 
permutations. In the usual approach one pick up a particular solution and the symmetrization 
is done explicitly when one computes the observables. For example one writes 

< Qa,b >= ry Z^5'n(a),n(fe); (1) 

where 11 denotes a permutation, 11 (a) is value of the index a after the permutation and the 
sum is done over all the n! permutations of n elements (n — i> at the end). 

This symmetrization is crucial when we need to compute quantities in finite volume because 
in this situation replica symmetry must be restored. The physical interpretation of sponta- 
neously broken replica symmetry as the coexistence of infinite many pure states is deeply 
rooted on this procedure. For example in the case of spin glasses we have that 

— J2 < ^h^^.■■■^^, >'-< Qib >= -jE?n(a),n(6)- (2) 

^^ ii,i2,...,ip=^,N "'■ n 

Things however become more complicated in systems which have one direction much longer 
that the others. Let us consider a system with NL degrees of freedom (e.g. <Jk,t), where 
k = 1,...,N, t = 1,...,L) and let us suppose that the interaction has a finite range in the 
variable t. We are interested in the case where L is very large (infinite) and A^ is large, but 
definitely much smaller than L. General arguments tell us that the correlation length in the t 
direction must be finite as soon N is finite: the system becomes one dimensional and in this 
situation no spontaneous breaking of a symmetry is possibile. 

For large A^ we can formally associate an order parameter to each system at fixed t. Systems 
at quite different t will have order parameters which point into different directions in the order 
parameter space (in our case replica space); indeed standard arguments imply that also if we 
constrain the order parameter at time zero to point in a given direction the symmetry will be 
restored at large times. This restoration is related to tunneling events among configurations 
in which the order parameters have different orientations in replica space. In this situation the 
zero modes connected to the spontaneous breaking of the replica symmetry are lifted and a 
finite correlation length (in the t direction) appears in the problem. Of course this correlation 
length diverges in the limit when N goes to infinity. 

This kind of problems appears in many cases, for example if we study the dynamics of spin 
glasses and we are interested in the tunneling among different pure states. A similar problem 
arises also in the equilibrium behavior of samples that are much longer in one direction that 
the other. 

In this paper we study these problems in a different context, the equilibrium behavior of 
infinitely long directed polymer. In the replica approach one finds that at low temperature 
the replica symmetry is broken. This breaking is an artefact of the mean field approximation 



and replica symmetry is eventually restored in this problem. In this case it is crucial to find 
out the solution to the problem of tunneling among different ground states, which is likely to 
control the large time behavior of the correlations. 

The number of degrees of freedom of a directed polymer embedded in a D- dimensional 
space is proportional to D. The replica approach can be successfully used when D goes to 
infinity, but for fixed D we have to deal with the existence of tunneling effects and replica 
symmetry restoration at large times. It is quite possible that this effects are responsible for 
the inaccurate value of the exponent computed at low dimensions in the replica approach. 

In this note we present a tentative computation of tunneling effects in the case of breaking of 
the replica symmetry in a continuous way, i.e. the function q{u) is continuous [|. We will show 
that one can define infinitesimal permutations and the tunneling among configurations differing 
by an infinitesimal permutation is not exponentially suppressed and gives a contribution that 
decreases as a power of D at large dimensions. 

In the general case the computations are rather involved. In the present case the analysis 
is simple because we can have to find the ground state of the Schroedinger equation for n 
interacting particles. The effects of the tunneling among wave functions with different order 
parameter may be taken into account in a first approximation by using a symmetrized wave 
function and computing the effects due to the superposition of different wave functions. 

This paper is organized as follows: in section 2 we recall some of the results of the replica 
approach to directed polymer, we present in the next section the symmetrized trial wavefunc- 
tion and we compute the corresponding free energy. The section 4 is devoted to the discussion 
of the integral over permutations and its consequences to the finite-D corrections. The last 
section summarizes the results obtained. 

2 Mean-field solution 



The directed polymer in a random external field is described by the Hamiltonian H] [ttO] [11 



%] = /dt(^(^) +V{tMm (3) 

where t is the coordinate along the polymer and the D-component vector uj describes the 
transversal coordinate. The random external field is supposed to be Gaussian-distributed with 
the correlation 



r(t, Uj)V{t', CO') = -DS{t - t')f{^-y=^). (4) 



The relevant part of the correlation is the long-distance behavior, which we choose to be 
power-law 

^^"^^ ^ 2(1^"^'^'^"^ , x ^ oo . (5) 

In this approach the directed polymer is a classic-mechanical static object, so that the 
interesting properties are related to its shape. Moreover, we neglect the size of the particular 
monomers. In this way there is no intrinsic natural length scale in the system. 

The scaling exponents are crucial to describe the large distance behaviour of the system. 
In particular ( controls the transversal fiuctuations 

{{uj{t) - uit'W) ~ (t - tr'^. (6) 



^In a separate publication we will discuss in more details the simpler case where one the replica symmetry 
is broken at one step, i.e. the function q{u) takes only two values ||8|. 



Sometimes ( is called wandering exponent. 

The exponent x controls the sample-to-sample fluctuations of the free energy of a finite 
part of the polymer, of the length L 



F(Ly-F{L) ~L'^. (7) 

There is a scaling formula relating ( and x? resulting from the Galilean invariance: x = 
2C-HT3 0. 

Different approaches to the calculation of the scaling exponents are possible [ill] Ell 



1^], among them the replica method is one of the most powerful. Making use of the replica 
trick, the problem is equivalent to the system of n particles, where n is the number of replicas, 
with pairwise interaction determined by the correlation function of the random external field. 
The exact solution is known for two opposite cases. For D = 1 and (5-correlated disorder the 
Bethe Ansatz leads to the value C = f- 

On the other hand, simple scaling argument (see e.g. [jl6|) leads to the Flory result 



2(1 + 7)' 

The same scaling analysis gives us the relation for the effective value of 7 for the case of 
(5-correlated disorder 7 = 1 -|- D/2, thus the Flory result does not coincide with the exact 
solution for D = 1. The question of the validity of the Flory formula arises and it is partially 
answered by the results obtained in another case, where the exact solution is known, namely 



in the D = 00 limit |Tl]]. In this case, the mean field approach should give the correct result. 
The mean-field equations can be obtained in a number of ways. In the most transparent 
method one consider the perturbation expansion for the propagator Gab{t — t') = {^ait)^bit')) , 
where a and b are the replica indices. The role of the order parameter is played by the self- 
energy aab defined by 

Gabit) = {[-§^ - cr]-\, (9) 

The point is that the Hartree-Fock approximation as exact in the D = 00 limit {D corre- 
sponds to the number of field components), thus giving the desired mean-field equations for 
a. Only cactus diagrams survive in this limit, thus reproducing the Hartree-Fock expansion. 
Non-cactus diagrams give formally corrections which are proportional to 1/D. 

Two cases are possible. 

• In the situation called 'short range correlations', the matrix a is either replica-symmetric 
or it corresponds to one step of replica-symmetry breaking. This solution is thermody- 
namically stable for 7 > 2. The resulting wandering exponent is C = |, i e. no influence 
of the disorder is seen. 

• For 7 < 2 another solution should be looked for and the natural choice is to break the 
replica symmetry in the hierarchical manner already familiar from the spin-glass theory 
[|l[. The matrix aab is then parametrized by the pair {a,a{u)), where a{u) is a function 
on the interval (0,1) and the number a is the diagonal element of the matrix aab- At the 
end we arrive at a surprising result stating that the wandering exponent has exactly the 
Flory value. 

The problem can be formulated in an alternative way by using the Feynman-Katz represen- 
tation for the solution of the Schroedinger equation. In this case, the usual replica approach 
is equivalent to finding the ground state of the imaginary-time quantum Hamiltonian. 



If we use the variational principle and we suppose that the trial wavefunction {\ip >) is 
Gaussian, i.e. 

{uj\tlj) = exp{-Q-,^uJaUJb), (10) 

we recover the results of the previously described replica approach. This approximation cor- 
responds to taking into account only one of the minima of the free energy and supposing that 
this valley has parabolic shape, i.e. it corresponds to the D-dimensional harmonic oscillator, 
whose ground state wavefunction is Gaussian. The matrix Q~^ of variational parameters is 
then equivalent to the self-energy in the Hartree-Fock approach. 

3 Symmetrized wavefunction 

3.1 General considerations 

In case where the transverse spatial dimensionality D is finite, the trial wavefunction should be 
generalized in such a way that all the equivalent valleys are taken into account. The situation 
is analog to the double-well potential where the barrier between the valleys are proportional 
to D. For D finite we have to include tunneling from one valley to the other one. The relation 
between the equivalent valleys is provided by the symmetry of the problem, in our case by the 
permutations of the n replicas. We denote Sn the group of all such permutations and S^ the 
subset of Sn obtained by excluding the identity. 

The problem we want to solve is to find (in the limit n —* 0) the ground state of the 

Hamiltonian 

1 " (9 1 

^ = -oE(7r-)' + ^E/(^K-^^l)' (11) 



2 1 duj„ ^. J~D 

where the attractive interparticle potential can have two forms 

f{x) = -gSix) (12) 

or 

The scaling analysis suggests the first case corresponds to the second case for 7 = 1 + D/2. 
We will use the form ([l3|) hereafter. 

In the usual approach one minimize the free energy functional 

F- 1 (^W) (14) 

nD {iPlij) ^ ' 

in the space of Gaussian wave functions. Here we want to generalize this approach by min- 
imizing the same functional in the space of the symmetrized wavefunctions of the Gaussian 
type 

\^) = y: niV'i) (15) 

nes„ 

{u\iJi) = exp(-- J2 Qab^a ■ UJb) (16) 

a,b 

with the matrix Q as a variational parameter. This choice corresponds to the prescription used 
in [|ll[] , where the role of Q~^ plays the selfenergy, which is constant in longitudinal momentum 
space. The permutation 11 of the replicas acts as a transformation of the matrices Q and Q~^: 

Qab ~^'^Qab = Qu{a)U{b),Qab ^^Qab = Qn(a)n{b)- (1^) 



Unfortunately enough when we try to compute quantities hke {'il)\H\%l)) or even {ip\ip)i we 
have not been able to write in a closed form the sum over the permutation in the case where the 
matrix Q is such to break explicitly the replica symmetry. It is quite possible that a compact 
formula for such a sum does exists and can be found in a smart way. 

In this note we have esplored the possibility of selecting those permuations we believe are 
the most relevant are we have done an approximate evaluation of the quantity defined in eq. 
p. 

3.2 Infinitesimal permutations 

Let us denote 5Q^^ = ^{IIQ^^ — Q^^) the difference between the permuted and the original 
matrix. Naively speaking we expect that for most of the permutations, for which 6Q~^ is a 
quantity of order 1, the contribution to eq. (|l^) is very small for large D, because the overlap 
between the original and permuted wavefunction tends to zero as the transversal dimensionality 
D goes to infinity. This is certainly true as long as the number of replicas remains greater than 
one. 

There are technical difficulties arise when we perform the replica limit n -^ 0. One of the 
peculiar things happening here is that the Schwartz inequality is no more valid (the space of 
the Q has negative dimensionality when n < 1) and in many cases it gets reversed. Sometimes 
one finds that overlap of the original and permuted vector blows up instead of vanishing when 
D ^ oo. However, we can relate this problem to the set of strange phenomena in the replica 
method such as the change of the mimima of the free energy into maxima which arise as an 
effect of the analytic continuation in n. 

The common procedure is to take the n ^ limit after all other calculations, which 
involves some changes of order of limits. In this case we will make the Z? — >■ oo limit first and 
the replica limit {n -^ 0) after that. More specifically, we will compute the integrals of the 
type J dxexp{a{n)Df{x)) by saddle point method in the region n > 1, where a is negative 
and then we will continue the result to ra < 1, where a is positive and strictly speaking the 
integral is not defined. However, in fact, the integral represents the sum over negative number 
of terms, thus the 'undefined' integral should be actually defined by the procedure we use to 
compute it or by another methodologically similar procedure. 

If we exclude large permutations, we expect that only the infinitesimal permutations, i.e 
those for which SQ~^ is small, will contribute. These infinitesimal permutation are not present 
for integer n and they arise only as an effect of the analytic continuation in n and of the 
peculiar way of breaking the replica symmetry. 



The infinitesimal permutations we will sum over were already classified by Goltsev |T^ . In 
his classification the exchanges of two replicas play crucial role as generators of the whole group 
of permutations. Clearly, the pairwise exchanges are not the only infinitesimal permutations; 
one should include the permutations of three, four etc. replicas as well, because we do not refer 
to the set of generators but to the set of permutations with small 6Q~^ in general. Nevertheless, 
we will use the pairwise exchanges of blocks of replicas as typical representatives of the whole 
set of infinitesimal permutations with the hope that they cover that set densely enough to give 
reasonable results. 

We suppose that the matrix Q~^ has the usual hierarchical form [|1^ 

Qaa=Q = (lW (^o^ 

Q-,' = q{x),xeiO,l) ^'^^ 

where x = a fl 6 is the overlap (ultrametric co-distance) between the different replicas a and b. 



Thus even the matrix Q, which is the inverse of Q ^, has the same hierarchical structure 

Qaa = r = r{l) , . 

and exphcit relations between the pairs (g, q{x)) and (f, r(x)) can be found in the appendix of 
[Plj . We will frequently use these relations in the last part of the calculations. We use double 
notation for the value of the diagonal element for convenience in the following calculations. 

Let us now define the permutations we will take into account. We will take two blocks of 
replicas of the size tub and of the mutual co-distance m. Clearly m < ms- The permutation 
in question, which we will denote 112 (m^, f^^) will involve the exchange of these blocks as fixed 
units. We will use the notation S':^'^ for the set of all permutations of this type, the identity 
excluded. 

Denoting the chosen blocks by the indices a and (3,a ^ (3 and the single replicas within 
these blocks by pairs (a, i) and (/3, j) it means that the following relations hold for the overlaps 



(a,i) n (/3,j) =mWi,j 

(a,i) n (a,j) = {/3,i) (1 {/3,j) = tub Vz,j 

The explicit expression for the change in the matrix Q~^ is 

+Sal3Sba^bl3{Q7ai)(bi) " Q70i)(bi)) + 



(20) 



+0ba0aaSal3{Q (^ai){l3j) Q{ai){aj))~^ 
+Sbf3SaaSal3{Q7ai)(ai) ' QjaiMBj))) ■ 



(21) 



We will rewrite the expression for the free energy 

nD (^i|^i) + Ene5i(V'i|V'n) 



f22) 



where \ipYi) = niT/^i) and \ipi) was defined in (|TB|). We now face the problem of evaluating 
the various terms in the previous equation as function of the infinitesimal permutation and 
eventually to sum over all the infinitesimal permutations. 

3.3 Some explicit formulae 



Here we evaluate the various terms which appear in eq. (|2^ and we postpone the sum over 
the infinitesimal perturbation to the next section. 

For the overlap between the original and the perturbed wavefunction we find 

(V^il^n) = det -^/2(g-i(i ^ Q5Q-')) = (V^i|^i) exp(-^rn) (23) 

where the function Tn can be expanded in powers of 6Q~^. 

Tn = TrQSQ-' - ^TtQ5Q-'Q5Q-' + ... (24) 

This power expansion will be the starting point of the following computations. Of course, 
the power of 6Q~^ does not tell us what is the order of smallness of the certain term, because 
the magnitude of SQ~^ itself depends on the parameters of the permutation in question. Thus, 
we will proceed by choosing an infinitesimal permutation and then collecting the terms of 



the same order of magnitude, which will arise during the computation of the value of the 



subsequent terms in (|2^) 



Just the same consideration holds for the matrix elements of the Hamiltonian between the 
original and permuted state. The kinetic part of the Hamiltonian gives 



(^n|-|EaAa|^l) = 

f [Trg-1 + \Tr6Q-'Q6Q-' - \Tr6Q-'Q6Q-'Q6Q-' + . . .]{^n\^i) 
and the interaction term 



(25) 



JZziQaa + Qbb — 2Q 



ab 



|l-7. 



{Qaa + Qbb - 2Qab)-H[QSQ-'Q]aa + [Q6Q-'Q]tb - 2[Q5Q-'QU)+ (26) 

{Qaa + Qbb - 2Qab)-mQ{6Q-'Q)Xa + [QiSQ-'Q^U- 

2[Q{6Q-'Q)Xb) + ...}{^Pu\^Pi) 
The first terms in each of these expansions correspond to the mean field {D = oo) approach 



IS 



Ho = Fd=oo = \q+ \9t^ f dx(2(f - r{x))Y~^ . (27) 

2 2 1 — 7 JO 

As we will see in the following, the lowest order of magnitude in the smallness of the per- 
mutation corresponds to the first and second order in 5Q~^, i.e. the terms up to the order 
0{{5Q~^Y) are necessary. 

The rescaled interaction constant g arises in the calculation. It depends on the bare coupling 
g and the transverse dimensionality D. For D ^ oo the two coupling constants coincide, g ^ g 



n]]. Because the interesting effects do not regard the rescaling of the interaction constant, we 
chose g independent of D and g having the corresponding D-dependence. 

When computing the replica summations involved in (^3p and (^) we use the tree diagrams 
which describe the configuration of replicas imposed by the ultrametric structure of the replica 
space. Let us turn to the value of Tn first. The lowest order in 5Q^^ corresponds to two 
non-zero trees shown on the figure la. The corresponding value is then 



TTIB 



TrQ6Q = 2mB / dx{r{x) — r{m)){q{x) — q{m)) . (28) 



Similarly, the trees contributing to the following term in the expansion of Tn are depicted on 
the figure lb,c and the value is 

-\Tr{Q5Q-^f = -{mB C dx{r{x) - r{m)){q{x) - q{m)f+ 

^sir - J^g r{y)dy - mBr{Tn)) 

[Im^ dx(f — J^ r{y)dy — xr{x)){q{x) — q{m)Y— 

C dx/™^ dy{r{mm{x,y)) - r{m)){q{x) - q{m)){q{y) - g(m))] . 



(29) 



We will suppose that the function q{x) has the same form as the one found in the case 
D = oo, i.e. 

q{x) = Ax^ X G (0,Xc) 

q{x) = Axl X G (xc, 1) (30) 

Q' = /^ + /o q{.x)dx . 

The parameter /i was introduced in order to have all quantities finite. It corresponds to 
enclosing the system to a finite volume of the linear dimension l//i. Using the relations 



between q{x) and r{x) listed in |TT| we find for m^ < Xc 



1 jY) Tfi 

Tn = T{ms, m) = --{s + 1)^(1 f + 0((1 )') (31) 

3 TUB niB 

while for itlb > Xc 

Tn = -i^ - 1 1 - - ' + 1 - -f . 32 

C C C C 

Hence, we see that the values of m which dominate the sum over the permutations are those 
near (and smaller than) itib for tub < Xc and those near Xc for tub > Xc. If we suppose 
the matrix elements of the Hamiltonian to have the same power dependence on 1 — m/mB 
and on 1 — m/xc in the respective regions of ttib, the dominating D-dependence of the free 
energy comes from the permutations with rriB < Xc- In the following we will see that the 
matrix elements of the Hamiltonian actually do have this property. Thus, the sum over the 
permutations means an integral over the parameters m, rriB, mB G (0, Xc) and m G (0, m^). 

Similarly we compute the terms arising in (ipulHltpi) . Summing of all tree diagrams gives 
for the kinetic part 

\Tr{6Q-yQ = H,{mB,m) = 



^{r - J^^ r{y)dy - mBr{m)) /™^ dx(g(x) - g(m))2 
+ ^[C dx(/J r(y)dy + xr(x) - f)(g(x) - g(m))2+ 
/"^dx/™^d?/(r(min(x,y)) -r(m))(g(x) - q{rn)){q{y) - q{m)) 



(33) 



Using the form of q{x) defined in (|30D we obtain, for m^ < x, 



H.irriB, m) = j-As{s + l)2m%+\l - —f + 0((1 - —f) . (34) 

48 TUB fUB 

For the interaction terms, similarly 
H2{mB,m) = 

YaHaAQ-a + Qbb - ^QabY-^msQ-^QU + [Q5Q-^QU - 2[Q5Q-'QU) = 

^2^"^{/^ dx /^ dy[(e(x, x) - 0(x, y)){r - r(min(x, y))Y^+ 

(9(x, x) + 6(x, y)){f — r(m))~'^] + 

!l dx[(20(x, 1) - 0(x, x) - 6(1, l))(f - r(x))-T- 

(2e(x, 1) + e(x, x) - 6(1, l))(f - r{m)Y^] + 26(1, l)(f - r{m)Y^} 



(35) 



where we have used the auxihary quantities 

e(x, y) = ^{%{x, y) + $o(y, x) + $i(a;, y) + <l>i(y, x)) 

$0(2^; y) = (^(a;) — r{rn)) J^^dz {q{x) — g(m)(r (min(z, y)) — r(m)) 

(36) 
$i(x,i/) = 

= (r(a;) — r{m)){q{y) — q{m)){yr{y) + Jydzr{z) — f) for y G {ra.mB) 

= otherwise. 

In the last term we will write only the lowest order (in {itlb — rn) ) terms. 

-\9 Ea^Qaa + Qbb - 2Qabr^ {[Q{SQ~'Q)Xa + [Q{SQ~'Q)\b " 2[Q{6Q-'Q) 



ab 



2 

-g2^-^^{J^ dx J^ dy [{{r{x) - r{m)) 



2 



{r{x) — r{m)){r{y) — r{m))){r — r(min(x, y)))~'^+ 

{{r{x) — r{m)y + {r{x) — r{in)){r{y) — r{m))){f — rm)^'^] + 2(f — r{m)Y^"' + 
Smdx [{f — r[x))^'^'{2{r[x) — r[m)){f — r{m)) — {r[x) — r[m)Y — (f — r(m))^)- 
(f — r{m))~"'{2{r{x) — r{m)){r — r{m)) + {r{x) — r(m))^ — (f — r(r«))^)]} 
/™^ dx (xr(x) + /^ dyr{y) — r){q{x) — q{m)Y + 0{{mB — mY) . 
Substituting the functions q{x) and r{x) we get the following results 

H^imsM = ^2-n-i^)-n2igi§y[(;^ - ^ + f±f ^)- 



-B 



Hs{mB,m) = -g2^-'r^{J^^ dx[^\x) /^„ d#-^(y) + (1 - a:) 0^ (a;) 0-7 (a;) - 



(37) 



^b{^{1 - gm^) - ^(1 - f±fx,))] + ^(1 - 7 + mB);;^}x (3J 

(dr. - ;pW + f±f ^)-^|(5 + 1)^(1 - ^r + 0((1 - m/m^)^) 



<j>{x) /^^ d#(i/)r^(y) - 0(a;)r^(a:) /^^ d#(i/)] + 

rn^i?) /m, dx[(l - mB)(P\x) + 0(x) /^^ dycf^iy)] + 202-7(^^)+ (39) 

/^, dxHrU^) - rH^B))i2<Pix)4>{mB) - 4>\rnB))- 

(0-^(0:) + r"(mB))02(a;)]}ls(s + l)m|5i(l - ^f + 0((1 - m/rriBr) 



where 

(j){x) = r{x) - r{mB) = —^^{^ - ^) 

(f)(x) =r- r{x) = - A'sis+2)(lF^ ~ ^ + ^^ 



(40) 



4 Summing over permutations 

Having computed all the matrix elements, we could proceed to computation of the free energy 



(|22D . For this end we first need to know the measure d/in in the space of permutations. Starting 



with the discrete p-adic formulation of the replica symmetry breaking [^], where the size of 
the i-th block is rrii = p* and p^ = n, we see that there are 

AP{mB, m, = ^ \ ^ ' 41 

permutations which exchange blocks of the size ms at the co-distance rrii. The continuum 
limit corresponds to taking p — > 1~, while K ^ oo and n = p^ -^ 0. Note that the usual 
replica limit n — > and the limit which introduces the continuous replica symmetry breaking, 
p — »• 1~ should be taken independently, as in the spin-glass theory |jl]]. Thus, the measure in 
the set S^^ of the permutations in question is 

d/i(^) m^,m =- ^, 3 42 

ImpvaB 

and it is immediately seen to diverge for p — > 1 ~ . 

The meaning of the divergence clarifies itself once we take into account the rest of the set S\ 
of all possible non-identical permutations. In the limit n— >-0|S'^|=n! — 1 = and the infinite 
number of permutations G S)f^ should be canceled by the contribution of other permutations. 
Actually, if we find e.g. the measure corresponding to the space S^^^ of all cyclic exchanges of 
three blocks of replicas, we see that their number is infinite again, but with opposite sign than 
the number of permutations in S^^ . It is that cancellation of infinities that makes finally the 
total number of non-identical permutations to be zero. 

It would be very difficult task to perform the same calculations as those presented above 
for all the variety of possible permutations and to show the cancellation of the divergences 
explicitly. Instead, we will approximate the integral over all permutations with the exact 
measure by the much more simple integral over the set S"^^^ only, but taking a "renormalized" 
measure dfi^'^\mB,Tn) instead of the true dix^'^\mB,m). The interpretation we attribute to 
such a replacement is that we pick out the permutation 11 with given ttib and m and suppose 
that the other permutations fl, exchanging not only two, but three, four etc. blocks of replicas 
are close to 11 in the sense that corresponding Tn and the matrix elements of the Hamiltonian 
are approximately equal, but when summing over all ll's, the divergences cancel, yielding the 
average measure dfi^'^\mB,rn) which is free of the divergence at p — >■ 1. Clearly, the above 
supposition is somewhat unfounded and its mathematical formulation is worth of further study. 
Nevertheless we consider the hypothesis to be plausible and useful at least for the purpose of 
establishing the power of the finite size corrections, which is what we are trying here. Thus, 
we will rely on it in the following. 

The question then arises, how to choose the approximate measure dfl^^\mB,m). The 
measure has to meet the condition 

/ dfx^^\mB,m) = dfj.^^\mB) (43) 



where dfi^^\mB) is the number of all permutations of block of the size between tub and 
rriB + drriB with the restriction that they cannot be expressed as a permutation with another 
(smaller) block size. This restriction avoids multiple counting of the permutations. Clearly, 
the sum over the block sizes from 1 up to given rriB should give {n/mB)\ , so 

d^(i)(mB) = T^r(— + l)dmB "=° -T'il)ArdmB . (44) 

dniB fTiB rriB 



Similarly, the following relation holds for the original measure /x 



(2) 



n 



dfi^'^\mB,m) = -— ^dniB (45) 

mG(o,ms) ZinprriB 



hence, per analogiam we infer that 

dpiP''{rnB-,'rn) = —V'{1) — ^dmdms . (46) 



The same result comes out also following another consideration. Dividing all the permu- 
tations into disjunct groups according their characteristics itlb (by the same means as we did 
when we were deriving /i*-^-*) and within a selected group taking only such a permutations 
which exchange two blocks, we obtain exactly the same result. This approach is based on 
the tacit assumption, that the pair exchanges dominate all quantities of interest. Actually, 
when we perform the calculation of Tn = T{mB-,Tni,m2) in the case of a cyclic exchange 
of three blocks (the permutation is characterized by two co-distances mi < 1712) we find that 
T{mB, TTLi, 7712)— T{mB, itli) ~ {l — {m2/'mB))^ > and for large D the pair exchanges actually 
dominate. 

The sentence just stated is slightly incorrect in that it holds only for n > 1. In the replica 
limit n -^ the inequality reverses and one would naively expect that, on the contrary, the 
pair exchanges are negligible. But the same consideration as the one which has clarified the 
wrong sign in the exponential exp(— DTn/2) applies here and the previous statement about 
the dominating permutations is restored. 

Finally, we obtain for the free energy 

F = Ho-r{l) f ^^!!^^H{mB,m)exp{-^T{mB,m)) (47) 
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where H = Hi + H2 + H-^. Note that the kinetic term. Hi is of higher order in {itib — itl) than 
the interaction one, H2 + H3, so the kinetic term can be omitted completely in the computation 
of the finite-dimensionality corrections. 

Setting rh = 1 — m/mB we can see that the free energy has the form 

F-Ho= f dmBCimB) f dm m^e^^^- 3)(^+i)'™' (48) 

and the integral over fh can be easily computed, yielding 

/' dmmV^('3)(-+i)'™^^ ~ _35/324/3r(l)(^ ^ IY^I^D~*'\ (D ^ 00); . (49) 
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The main inference from this result is the D-dependence of the correction to the free energy, 
which is 

F-Hor^ D-^'^. (50) 



Of course, the full solution should continue with minimizing the free energy with respect of the 
matrix Q^^, but the non-integer power of D in the corrections to free energy and propagator 
must persist. The open question and the most interesting one is, however, what are the finite-D 
corrections to the exponent s, which is the quantity which governs the long-distance behavior 
of a directed polymer. In answering this question, we cannot avoid the explicit computation of 
the stationary point of the free energy. In principle it is straightforward: we take the formulae 
just obtained and let the first derivatives with respect to the three variational parameters 
s, A, Xc be zero. This work is in progress and we will refer the reader to further publication. 

5 Conclusions 

We have seen that infinitesimal permutations can be identified in the case of continuous break- 
ing of the replica symmetry. The sum over these infinitesimal permutations is not a simple 
task: if one proceeds in a naive way the measure over a class of infinitesimal permutations is 
formally infinite and this infinity is related to having limited our attention on a restricted class 
of permutations. The corresponding double counting problem gives rise to this divergence. We 
tentatively propose to cure this divergence by introducing an effective measure in order to take 
care of these double counting problems. 

If we follow this tentative approach we can compute the shift in the ground state energy 
due to the presence of the tunneling and we find that in a given range of the parameters there 
are corrections which go to zero like D~^'^. The effects of these corrections on the critical 
exponents has not been studied yet. However, we think that the tunneling among different 
states plays a crucial role in the final theory. It would be extremely interesting to compare these 
results with those coming from a standard 1/D expansion to see if these effects are correctly 
taken account in the perturbative 1/D expansion or they have a more non perturbative nature, 
as this study may suggest. 
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Fig. 1 a, b and c. Tree diagrams corresponding to the summations over 
rephca indices in the equation (|2^). 



